Abstract. We present a fit to the spectral energy distri- 
bution of OH 127.8+0.0, a typical asymptotic giant branch 
star with an optically thick circumstellar dust shell. The 
fit to the dust spectrum is achieved using non-spherical 
grains consisting of metallic iron, amorphous and crys- 
talline silicates and water ice. Previous similar attempts 
have not resulted in a satisfactory fit to the observed spec- 
tral energy distributions, mainly because of an apparent 
lack of opacity in the 3-8 fj,m region of the spectrum. Non- 
spherical metallic iron grains provide an identification for 
^vq ] the missing source of opacity in the near-infrared. Using 

■ the derived dust composition, we have calculated spectra 
' for a range of mass-loss rates in order to perform a consis- 
• tency check by comparison with other evolved stars. The 

P5 L — [12 fim] colours of these models correctly predict the 
^ ^ ■ mass-loss rate of a sample of AGB stars, strengthening 
[ our conclusion that the metallic iron grains dominate the 

■ near-infrared flux. We discuss a formation mechanism for 
I ' non-spherical metallic iron grains. 
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1. Introduction 

Oxygen-rich Asymptotic Giant Branch (AGB) stars show 
an infrared excess on their spectral energy distribution 
(SED), which arises from thermal emission of dust located 
in a circumstellar shell. It is well established that the dust 
in this shell mainly consists of silicates, deduced from the 
clear detection of the 10 and 18 /xm features. These bands 
are due to the Si-0 bond stretching and Si-O-Si bond 
bending vibrations, respectively. It was noticed already 
some 



Le Sidaner & Le Bertre 1993). Recently there have been ef- 
forts to derive improved optical constants for astronomical 
silicates, both for lines-of-sight in the interstellar medium 
(ISM) and toward oxygen-rich evolved stars exclusively 
( psscnkopf et al. 1992| ; [David fc Pcgouric 1995| ^uh 1999| ) . 
These optical constants still do not provide the answer to 
the composition and the nature of circumstellar silicates. 

The study of circumstellar dust shells in the pre-ISO 
era was based on low resolution spectroscopy, but with the 
launch of t he Infrared Space Observatory (ISO) ( Kesslei 



that the opacities of vufluus Lypes of slUcaLes measured lu 
the laboratory, were not high enough to explain the shape 
of the spectrum in the near-infrared (NIR, 3 < A < 8 fim). 
Jones &: Merrill (1976| ) investigated the effect of pure sil- 
icates with graphite inclusions and of meteoritic rock on 
the shape of the SED. These materials were taken to repre- 



^5 years ago JlJoiics fc MerrUl 1970|; [Bedljii 1977]) ^^^^ ^^ ^ ^^^^^ wavelength region, from 2-200 ^m. 



sent siliicaLes wiLli 



wiLli various metalhe luelusiuus, and referred 
to as dirty silicates. Radiative transfer ealeulatiuiis using 
the opacity of dirt y silicates were in goo d agreement with 
the observations. Draine fc Lee (1984 ) used laboratory 
measu rements of limited resol ution and wavelength cov- 
erage (Huffman fc Stapp 1973) in which the optical con- 



became available for intermediate resolution spectroscopic 
observations (A/AA > 400). The ISO spectra of oxygen- 
rich evolved stars turned out to be extremely rich in 
solid state features. Hence, it became meaningful to at- 
tempt spectral fits using properly characterized cosmic 
dust analogs, rather than astronomical silicates. Demyk 



et al. (2000) aim to do this for OH/IR stars (i.e. AGB 



stars with an optically thick dust shell), but their work is 
not complete and parts of the spectrum are not fitted to 



stants are partly modified and expanded in wavelength 
regime to match the astronomical observations of inter- 
stellar and circumstellar silicates (Jones fc Merrill 1976; 



Rogers et al. 1983) These optical constants are usually 
referred to as astronomical silicate. The exact chemical 
composition of this astronomical silicate is not known, 
nevertheless it is widely used in theoretical studies of the 
radiative transfer in circumstellar dust shells (e.g. BedijE| 
19871; ^chutte fc Tielens 19891; Ljusttanont fc Tielens 19921; 
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a satisfactory level. Harwit et al. (2001) have studied the 
optically thin dust shell surrounding red supergiant VY 
CMa and are able to fit the SED, albeit after including a 
large overabundance of metallic iron in the dust composi- 
tion. 

In this study, we aim to determine the dust composi- 
tion with full radiative transfer calculations using exclu- 
sively optical constants of well-defined materials measured 
in the laboratory, examining all astronomically relevant 
materials, and taking abundance constraints into account. 
Revealing the exact dust composition is an important step 
toward understanding dust formation and processing in 
the outflows of oxygen-rich evolved stars. 

The paper is organized as follows: In Sect. || we dis- 
cuss the spectral energy distribution arising from full ra- 
diative transfer calculations using amorphous olivine as 
the only dust component, and compare it to the obser- 
vations. Sect. H discusses the dust composition and grain 
properties required to improve the fit to the observations. 
A spectral fit to OH 127.8-1-0.0 is presented in Sect. H, 
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supported by a consistency check on other AGB stars in 
Sect. |5[ Section |^ contains a discussion of the results. 



2. Modelling the circumstellar environment 

In this section we describe and discuss the results from 
the radiative transfer calculations assuming that the dust 
only consists of amorphous olivine. We describe the dis- 
crepancies between the calculated and observed spectral 
energy distributions of AGB stars. 



2.1. Model assumptions and default grid 

Our approach is to calculate a set of model spectra of 
dust shells characteristic for both Mira and OH/IR stars, 
and to compare these with ISO observations. The stellar 
and wind parameters are the same for all models in the 
grid, except for the mass-loss rate which is varied from 
relatively low values, typical for Mira's, to relatively large 
values, characteristic for OH/IR stars. This leads to a va- 
riety of optical depths, which can be compared to the op- 
tical depth in the 10 /im feature of observed spect ra. A 
similar approach has been taken by Bedijn (1987 ) who 



also studied the appearance of the SED of AGB stars as 
a function of mass loss, comparing the results with IRAS 
observations and ground based infrared data. 

We use the code MODUST to model the spectrum of 
the circumstellar dust shell. The radiative transfer tech- 



nique applied in this code has been outlined by Kemper 
(2001| ). The specification of grain properties, such 



et al. 



as size and shape distribution, are discussed in Bouw- 
man et al. (2000). Here, we only discuss details relevant 



for the presented models. The dust is assumed to be dis- 
tributed in a spherical shell, with inner radius i?in and 
outer radius i?out, which is irradiated by a central star 
with Teff = 2.7 • 10^ K, a radius i?* = 372 i?© and a lu- 
minosity = 6.3 • 10'^ Lq, thus describing a typical AGB 
star. The input spectrum of the central star is not simply 
a black body but is a characteristic mean spectrum for 
spectral type M9III, inferred from observations (380 nm 
< A < 900 nm) and extended with synthetic spectra in the 
range 99 - 12500 nm ( |Fluks et al. 1994| ). The M9III star 
with a Teff = 3126 K, emits significantly less flux at A < 1 
fxm than a black body of the same temperature. Dust is 
efficie ntly heated in this wavelength region, and therefore 



Table 1. Inner radius of the dust shell used in our initial 
model calculations, for different mass-loss rates, based on 
the work of pedijn (198^ ). 



M (Mq yr-M 


^in 


5(-8) 


6.05 


l(-7) 


6.20 


2(-7) 


6.37 


5(-7) 


6.56 


l(-6) 


6.91 


2(-6) 


7.26 


5(-6) 


8.23 


l(-5) 


9.17 


2(-5) 


10.2 


5(-5) 


12.3 


l(-4) 


14.3 


2(-4) 


17.8 


5(-4) 


23.5 


l(-3) 


29.0 



tion between the two possibilities, because only the dust 
in the outer regions is visible for the observer. 

The adopted mass-loss rates range between 5 • 10^* 
and 1 • 10^^ Mq yr^^. The inner radius of the dust shell 
is defined by the condensation temperature of amorphous 
silicate, taken to be Tcond = 900 K ( |Bedijn 1987| ). The 
exact location of the inner radius is determined by the 
largest grains and is a function of the mass-loss rate as the 
dust density in the shell determines the diffuse component 
of the radiation field, also known as backwarming, thus 
constraining the temperature profile. The values for i?in 
range from 6.05 to 29 i?* for the smallest and largest mass 
loss, respectively (see Tab. |l|). The outer radius is chosen 
to scale correspondingly to 100 R\a. This is significantly 
lower than what is used by others (e.g. Bedijn 1987), but 
higher values lead to an overestimation of the flux around 
20 - 25 /im. The dust density in the shell is determined 
using a constant outflow velocity of 20 km s~^, implying 
the density decreases as p{r) oc r^^. We used a dust/gas 
mass ratio of 0.01. 



To obtain the absorption and scattering efficien- 
cies we used optical constants o f an amorphous olivine 
(MgFeSi04) sample measured by Dorschner et al. (1995 ) 
and assumed that the grai ns are spherical, such that Mie 
theory can be applied (see van de Hiflst 1957 ; Bohren fc 



the dui ft located at the inner edge of the dust shell around 
an M9III star would be somewhat cooler than the dust 
at the inner radius around a black body of 3126 K. This 
temperature difference decreases with distance from the 
central star. In case of an optically thin dust shell, the ob- 
served spectrum is dominated by the thermal emission of 
the warm dust at the inner radius. Therefore, the silicate 
features are somewhat stronger in case the black body is 
used as a central source. In case of an optically thick dust 
shell, there is no difference in the spectral energy distribu- 



Huffman 1983 ). The size distribution of grains was as- 



sumed to be interstellar, i.e. given by n(a) oc , where 
a is the grain size ( Mathis et al. 1977 ). For the circumstel- 
lar silicates we adopted a minimum and maximum size of 
the grains of 0.01 and 1 /im. 

2.2. Problems in fitting the SED of AGB stars 

Figure ^ shows a comparison of two of our predicted 
SEDs with the ISO spectrum of the oxygen-rich star OH 
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Fig. 1. ISO spectrum of OH 127.8+0.0 (solid line), to- 
gether with two model spectra. The dashed line is the 
SED of the model with = 1 • lO"'^ Mq yr-\ and the 
dotted line corresponds to the model with M = 2 ■ 10"'' 
Mq yr"-'. Note the discrepancies between the model fits 
and the observed spectrum at A = 3-8 ^m region, at the 
peak position of the 18 /im feature and at the long wave- 
length part of the spectrum. 

127.8+0.0 (solid Hne). The models with mass loss 1 • lO"'* 
(dashed hne) and 2 • lO"'' Mq yr"! (dotted line) were se- 
lected. The first because it fits the 10 fiui region best; the 
second because it gives teh best overall shape of the spec- 
trum, even though the discrepancies remain large. The 
poor agreement between the model spectrum and obser- 
vations illustrates the problem of fitting AGB spectra. The 
main discrepancies are: i) the NIR flux is overestimated, 
ii) the shape and predicted position of the olivine reso- 
nance at 18 /Ltm are not well reproduced, complicating the 
fit to the mid-infrared (MIR, 8 < A < 30 ^m) spectrum. 
Specifically, we predict the center of this absorption fea- 
ture at 16.3 /xm. iii) The flux in the far-infrared (FIR, 
A > 30 /xm) is underestimated, and the models show a 
steeper spectral slope than observed. In Sect. |^ these prob- 
lems will be discussed and resolved. All three issues can 
be explained in terms of dust composition and grain prop- 
erties. 

3. Improvements to the spectral fit 

3.1. Metallic iron as a source of NIR opacity 

Let us first focus on the behaviour of the two predicted 



Fig. 2. Optical depth towards the central star in the dust 
shell of an evolved star with a gas mass-loss rate 1 ■ 10~^ 
Mq yr~^. The used optical constants are from amorphous 
olivine ( Dorschncr et al. 1995 ). r = 1 is indicated with 
the dotted line. The dashed line indicates absorption and 
the dashed-dotted line denotes scattering. The solid line 
represents the combined optical depth due to both effects. 
Note that the combined optical depth is dominated by 
scattering for A < 8 fim and by absorption for A > 8 /im. 



decreases with increasing density as a larger fraction of 
the thermal radiation emitted by the hottest grains suf- 
fers additional thermalization in colder regions of the dust 
medium. So, one way to further reduce the NIR flux is sim- 
ply to increase the mass-loss rate. However, this will also 
increase the optical depth at 10 and 18 ^m, and by that 
will overpredict the absorption strength of these features. 
Thus, varying the physical parameters, such as the den- 
sity distribution, the inner radius and the outer radius, 
will not cause a wavelength dependent effect on the opti- 
cal depth. In order to resolve the NIR problem one needs 
to selectively increase the opacity at these wavelengths, 
i.e. a species is missing of which the absorption efficiency 
follows a smooth distribution - as no unidentified reso- 
nances are present - predominantly contributing in the 
3 — 8 /im region. Using the model grid calculations to fit 
the optical depth at various wavelengths, we can estimate 
that TNiR/rio^m ~ 1/5. 

For chemistry and abundance reasons we considered 
the following dust species which could be present in 



oxygen-rich dust shells: metallic Fe (Ordal et al. 1988) 
FeO (iHenning et al. 1995|), Si02 ([Gray 1963|; prummond 



spectrE In both models the radial monochromatic absorp- 
tion oj tical depth in the NIR is about equal to or some- 



what larger than unity; in the 10 and 18 iim feature the 
dust medium is optically thick, while in the FIR it is opti- 
cally thin. This is shown in Fig. ^ where the optical depth 
is given for a typical dust shell consisting of only amor- 
phous olivine. In general, the flux at short wavelengths 



1936; Spitzer & Kleinman 1961; Philipp 1985; Longtin 
et al. 1988|), FC3O4 (|Stcyer 1974] ), AI2O3 ( [Koike et al. 1995| 



Begemann et al. 1997) and amorphous and crystalline 



H2O ice (jBertie et al. 1969; Warren 1984). The mass ab 



sorption coefficients and extinction efficiencies have been 
determined from the optical constants. The missing source 
of opacity has to contribute significantly to the continuum 
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in the 3-8 fj,m region, but should not show strong and nar- 
row features in the 2-20 region. This eliminates Si02 
and water ice for this purpose. Si02 has very strong res- 
onances at ^8.7, '^12.4, ^^14. 4 and ~20.5 /zm, while the 
contribution to the continuum in the 3-8 //m region is 
orders of magnitude less. These sharp resonances are not 
observed, therefore only a small amount of Si02 can be 
present in the dust shell, which is not enough to signif- 
icantly increase the NIR opacity. Water ice, both amor- 
phous and crystalline, has features at ~3.1, ~4.4, ~6.3 
and ~12 /xm. In between these features, the continuum is 
relatively weak. The ^^3.1 micron feature is the strongest 
feature and is actually seen in absorption in most AGB 
spectra. The amount of water ice required to fit this fea- 
ture explains only a negligible fraction of the NIR opacity. 
Fig. H shows the mass absorption coefficients of the re- 
maining candidates, with broad band or continuum emis- 
sion in the NIR wavelength region, calculated for both 
spherical and non-spherical particles. A continuous distri- 



bution of ellipsoids (CDE, see Bohren & Huffman 1983) is 
used to represent non-spherical grains. The strongest con- 
tribution in the 3-8 fim region is caused by non-spherical 
metallic iron particles, suggesting these are a likely candi- 
date for the missing opacity in the NIR in the spectrum of 
AGB stars. The mass absorption coefficients of FeO and 
AI2O3 are much lower in the NIR region, while the con- 
tribution at longer wavelengths is more significant. This 
rules these species out as the missing source of opacity. 
Fe304 does contribute in the NIR although a factor of 
2 less than Fe, but shows spectral features in the 15-20 
fim region which are not observed in the spectrum. In 
this study we will therefore concentrate on Fe as the most 
likely candidate for the moving NIR opacity, because it 
has the highest mass absorption coefficient in the 3-8 /im 
region and the most favourable spectroscopic properties. 

3.2. Non-spherical grains 

The shift in the peak position of the 18 /im feature between 
model calculations and observations is a second problem 
that occurs. This can be solved using non-spherical amor- 
phous olivine grains. The difference between the mass ab- 
sorption coefficients of spherical and non-sperical olivine 
grains is clearly shown in Fig. 2 of the work by Demyk 
et al. (20"00| ). In particular the shape of the 18 /im feature 
is much affected and the peak position shifts to a some- 
what longer wavelength for non-spherical grains. However, 
Demyk et al. (200C) did not use these non-spherical par- 
ticles in their radiative transfer calculations of optically 
thick dust shells. 

A continuous distribution of ellipsoids (CDE) may bet- 
ter approximate the actual variety of grain shapes that is 
present in the outflow of evolved stars. The optical proper- 
ties are dependent on the shape, which is for example ap- 
parent in a shift of the peak position of spectral features. 
Shape effects are most prominent for conducting mate- 
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Fig. 3. Mass absorption coefficients k (cm^ g~^) as a func- 
tion of wavelength for different species contributing in the 
NIR. The left panel shows the calculations for a continuous 
distribution of ellipsoids (CDE) which supposedly repre- 
sents non-spherical particles, and the right panel shows 
the Mie calculations, used for spherical particles. The la- 
bels indicate the mass absorption coefficient curves for the 
considered species. Note the large differences between the 
mass absorption coefficients for CDE and Mie calculations 
in case of metallic iron and Fe304, caused by their highly 
conductive nature. 



rials. Strictly speaking, the CDE approximation is only 
valid in the Rayleigh limit, where the size of the particles 
is at least a factor of ~ 20 smaller than the wavelength. 
As we are using CDE calculations to solve problems with 
the 18 /im feature and our maximum grain size used is 1.0 
/im, this condition is not violated. The result that metallic 
iron in CDE accounts for the missing NIR opacity, is also 
consistent with the Rayleigh limit, because the metallic 
iron grains are probably incorporated in the silicate grains 
(see Sect. ^) and therefore have sizes significantly smaller 
than 1.0 fim. For an elaborate description of resonances 
and the absorption characteristics of non-spherical parti- 
cles the reader is referred to Bohren fc Huffman (1983| ). 



3.3. Water ice features in the far-infrared 

The improved fit to the NIR part of the SED leads to a 
redistribution of the radiation such that the flux in the 
NIR is suppressed and FIR flux levels are increased. Ad- 
ditional improvement of the FIR flux might be achieved 
by adjusting the outer radius of the dust shell, where cold 
dust - optically thin in the FIR - contributes. The remain- 
ing discrepancy between the observed and modelled flux 
levels at the long wavelength side can probably be over- 
come by water ice as an additional dust component. The 
presence of the 3, 43 and 60 /im features in most OH/IR 
stars implies that crystalline water ice is an important 
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dust component ( Sylvester et al. 1999 ). The 3 feature 
is seen in absorption (see also Sect. 3.1). In addition to the 
43 and 60 /^m emission features, crystalline water ice also 
has an instrinsically strong underlying continuum > 30 
/xm which contributes to the FIR emission in the SED. 

4. Results for OH 127.8+0.0 

4-.1. fflpZ/Tf pnrnmp-fprQ nnrl rnrli/THiiP irnriQfpr mnfiplhrirj 



4.1.1. Dust properties 

We find that the dust shell consists of the following dust 
components: amorphous olivine (80% by mass), forsterite 
(3%), enstatite (3%), metallic iron (4%) and crystalline 
water ice (10%, Tcond = 150 K), based on the assumption 
that all dust species are present in separate grain pop- 
ulations. The crystalline silicates forsterite and enstatite 



are identified in the spectra of AGB stars (Sylvester et al 



1999), and their relative mass fraction is assumed to be 



We h ave constructed a model representative tor AGB 3% e ach, a typical value for AGB stars ( |Kempcr et al 



stars, [with the intention to fit the spectrum of OH 200l|) . We did not attempt to fit the actual mass frac 



127.8+0.0, a high mass-loss rate AGB star. The input 
spectrum at the inner radius is that of a star with spectral 
type M9III (|Fluks et al. 1994D , and a radius of 372 i?0 . The 



density distribution and total mass of the dust shell are 
determined by the velocity profile and the inner and outer 
radii. The outflow velocity was again set to 20 km s~^, 
implying a density distribution oc r~^. In reality the wind 
probably accelerates while dust is being formed in the in- 
ner regions, until it reaches the terminal velocity. However, 
the dust shell is very optically thick in the NIR, where the 
warm dust in the inner regions most effectively radiates. 



implyi] ig that density variations in the inner parts - well 
inside the r = 1 surface ~ of the dust shell do not affect the 
emerging SED. Therefore it is impossible to further con- 
strain the velocity profile and the inner radius; we chose 
the latter to be consistent with a condensation tempera- 
ture of - 900 K, using i?i„ = 14.3 i?* (3.7 • lO^"* cm). This 
temperature corresponds to the condensation temperature 
of amorphous silicates, but crystalline silicates and metal- 



tion taken by enstatite and forsterite, which is subject 
to further study (Kemper et al., in prep.). The overall 
shape of the SED is not affected by assuming a degree of 
crystallinity of 3%, since aside from the sharp resonances 
in the MIR, the opacities of amorphous and crystalline 
silicates are comparable. All dust is found to be present 
in non-spherical (CDE) grains. The adopted optical con- 
stants of crystalline and amorphous silicates are provided 
by Dorschner et al. (1995| ), Jager et al. (1998D and Koike 
(priv.comm.). For crystalline water ice we used the opti- 
cal constants derived by Bertie et al. (1969| ) and Warren 
(1984[); the optical constants of metallic iron arc taken 



lie iron condense at higher temperatures (Gail & Sedlmayr 



1999). We performed test calculations in which we varied 



the inner radius and did not see a difference in the re- 
sulting SED. It is possible to determine the outer radius, 
found to be at i?out = 5000 i?* (1.3 • 10^^ cm), with an 
accuracy of 15%, by fitting the slope in the FIR. By fit- 
ting the optical depth of the 10 /im silicate feature, the 
dust mass-loss rate was found to be 7(±1) • 1O~^M0 yr~^, 
which translates in a gas mass- loss rate of 7(±1) • 1O~^M0 
yr~^, assuming a dust/gas ratio 0.01. Together with i?out 
this parameter determines the total dust mass in the shell 
surrounding OH 127.8+0.0 and sets it to be 1.4-10"^ Mq. 

By scaling the SED to the observations, we find that 
the distance toward OH 127.8+0.0 is 1800 (± 50) pc, 
which assumes the stellar luminosity is 6300 Lq. This is 
significantly smaller than results using the phase lag of 
the 1612 MHz OH mascr in combination with imaging. 
Herman fc Habing (1985 ) find a distance of 6.98 kpc, and 
Bowers fc Johnston (1990 ) re port 6.21 ± 1.0 kpc, usin g 
the phase lag determination of Herman fc Habing (19851 ). 



Such a large distance consequently leads to a very high 
luminosity (2.6 • 10^ Lq), more typical for red supergiants 
than OH/IR stars. Other studies result in distances of 5.6 
kpc (lEngels et al. 1986|) and 2.8 kpc QHeske et al. 1990|). 



from lOrdal et al. (1988| ) . 

Fig. I shows the fit to OH 127.8+0.0. Th e ISO spec- 
troscopy is taken from Sylvester et al. (1999 ) and the 60 
and 100 /im photometry originates from the IRAS Point 
Source Catalog. Additional observations with the James 
Clerk Maxwell Telescope (JCMT) have been performed to 
obtain the 450 and 850 /im photometry points (Kemper et 
al., in prep.). It is clear that if 4% of the total dust mass 
is contained in non-spherical metallic iron grains, the NIR 
problem is solved. This is consistent with the average in- 
terstellar Fe and Si abundances ( Snow fc Witt 1996 ) that 
allow a mass fraction of 15% in the form of metallic iron, 
with respect to amorphous olivine. 

In order to fit the FIR part of the SED, it is required 
to include a significant fraction of water ice in the dust 
shell (10%) for which a condensation temperature of 150 
K is assumed. This is consistent with the mass fraction 
contained in water ice in the oxygen-rich dust shell of post- 
AGB star HD 161796 (Hoogzaad et al., in prep.). 

4-2. Comparison with astronomical silicate 

The derived total dust extinction toward OH 127.8+0.0 
can be compared with astronomical and laboratory mea- 
surements. The extinction efficiency Qoxt of the individual 
dust species assuming CDE have been added proportion- 
ately and the resulting total Qext is plotted in Fig. ^. Note 
that all curves are normalised to the strength of the 10 //m 
feature to emphasize the differences in the NIR region. To 
provide some quantitative insight. Tab. |^ presents the ex- 
tinction efficiencies at 7 /im and at 18 /im with respect 
to the extinction efficiencies at 10 /im. The long-dashed 
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Fig. 4. The large panel shows the fit (dashed line) to the ISO spectrum of OH 127.8+0.0 (solid line). In the inset, the 
same fit is presented for the 30 - 1500 fim wavelength range. The diamonds indicate the IRAS 60 and 100 /im and 
JCMT SCUBA 450 and 850 /im photometry points. The error bars on the measurements are smaller than the size 
of the symbols. A mass fraction of 4% of metallic iron was included in the dust. See Sect. |^ for a description of the 
model parameters and the dust composition. 



line shows the extinction efficiency for a synthesized pure 
amorphous olivine (Dorschner et al. 1995). One can see 
easily that for this species the NIR extinction is very small 
compared to the efficiency at 10 /im. Also shown in the 
figure ar e astronomical silicates derived from ISM lines- 
of-sight ( Drainc fc Lee 1984 ; Ossenkopf et al. 1992) and 



towar d late-type stars (pssenkopf et al. 1992; [David fc 



Pegoui ic 1995 ; 3uh 1999| ). Only the work of Osscnkopl 



et al. ( |1992| ) assumes non-spherical grains and uses CDE 
as a representation of grain shapes. The other efficiencies 
are calculated using Mie scattering for spherical particles. 



All in all, the David & Pcgourie result gives the poorest 
match to our findings. 3uh (1999) has derived the optical 
constants from AGB stars and splits them into two differ- 
ent types: warm and cool silicates, representative for low 
and high mass-loss rate AGB stars respectively. This di- 
chotomy assumes a difference in composition of the bulk 
material. In Fig. |[ we present the efficiency for cool dust 
as this is most appropriate for the optically thick dust 
shell modeled in this work. The Suh result compares rea- 
sonably well with our finding, though the NIR efficiency 



is less and the 18 /im efficiency is somewhat larger. Os 
senkopf et al. (1992| ) distinguish between interstellar 



The optical constants derived by David & Pegourie 
(1995|)]from G and M supergiants, closely resemble the 
synthetic olivine in the NIR region. The 10 /im resonance 
is much narrower, while the efficiency at longer wave- 
lengths is again significantly less compared to our result. 



rich silicates (olivines) and circumstellar 0-deficient sil- 
icates (pyroxenes), without discriminating between opti- 
cally thin and thick dust shells. Shown in Fig. || are the 
0-deficient silicates which represent the most meaningful 
comparison with our model. In the NIR the Ossenkopf et 
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Table 2. Relative extinction efficiencies at 7 and 18 /im compared to the extinction efhciency at 10 /im, in addition 
to Fig. I 



optical constants 
OH 127.8+0.0 

O-rich silicate 
O-deficient silicate 
astronomical silicate 
warm silicate 
cold silicate 
astronomical silicate 
amorphous olivine 



Qcxt(7^im) 



t(18Atm) 



Qcxt(10/xm) Qcxt(lOMm) 



ref. 



0.22 

0.15 
0.12 
0.087 
0.085 
0.080 
0.041 
0.017 



0.67 

0.61 
0.47 
0.39 
0.37 
0.78 
0.12 
0.68 



Dssenkopf et al. (1992 ) 



Dssenkopf et al. (19921) 



Draine fc Lee (1984) 



5uh (1999 ) 



5uh (1999D 



David 



(19951) 



Dorschner et al. (1995) 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 



0.0 t 



this work 
synthetic olivine 
Draine & Lee 
David & Pegourie 
Ossenkopf et al. 
Suh 




1 



1 5 



20 



X (/xm) 

Fig. 5. Comparison of the extinction efficiency Qext for 
the composition found for the circumstellar dust shell of 
OH 127.8+0.0 (this work), co mpared with the Qcxt o f syn- 
thesized amorphous olivine ( Dorschner et al. 1995), and 
the Qc xt for several different astronomical silicates (Draine 



fc Lee ^984 [Ossenkopf ct al. 1992| ; [David fc Pegourie 1995| ; 



Suh 1999| )?^A11 extinction efficiencies are normalised on the 
10 /im peak strength to allow comparison of the NIR ex- 
tinction. 

al. result compares best with our efficiencies, although the 
still find a lower Qcxt- In the 18 /im region their efficiencies 
also fall below our work. 

So, the dust composition determined in this paper pro- 
vides the only result that produces sufficient NIR flux. One 
should note, however, that our results are of course sensi- 
tively dependent on the metallic iron content of the AGB 
star, and may vary from source to source, whereas other 
studies tried to derive mean values for a sample of stars. 

5. A consistency check 

The SED of OH 127.8+0.0 is typical for high mass-loss 
rate AGB stars, which justifies the calculation of a series of 
models for varying mass-loss rates, while the dust compo- 



sition and other physica l parameters were kept the same. 
[Gail fc Scdlmayr (1999| ) show that metallic iron and sili- 
cates condense simultaneously in the stellar outflow, sug- 
gesting that the relative mass fraction of iron compared 
to the silicates remains constant. We allowed the inner ra- 
dius to change in order to satisfy the constraint that the 
condensation temperature Tcond = 900 K. 

Using the SEDs predicted by the model, the infrared 
colours in different stages of the mass-loss evolution on 
the AGB can be determined. For example, the relation 
between the K — [12 /xm] colour or K — L colour and the 
optical depth is a useful diagnostic tool to determine the 
mass- loss rate of an AGB star (Le Sidaner & Lc Bertrc 



1996; Le Bertre & Winters 1998). We searched for the 
presence of a relation between the L — [12 /im] colour in- 
dex and the mass-loss rate, since the L — [12 /im] colour 
can be determined from ISO SWS spectroscopy directly. 
For the L-band (3.6 /im) we use the transmission profile 
provided by Bessell & Brett (1988); for the magnitude 



calibration an abs olute flux density for zero magnitude of 
277 Jy was used ([Koornnccf 1983[) . The standard IRAS 
transmission profile and calibration are used to determine 
the magnitude in the 12 /im band. The calculated abso- 
lute magnitudes and the colour indices are summarized in 
Table |. 

Fig. H shows the correlation between the L — [12 /im] 
colour and the mass-loss rate derived from our model 
calculations. In addition, we have included the mass-loss 
rates and infrared colours of five stars taken from the sam- 
ple of Sylvester et al. (1999). The L — [12 /im] colour is 
obtained directly from the ISO spectroscopy. We can see 
that within the accuracy of the mass-loss rate determi- 
nations, a metallic iron fraction of 4% gives reasonable 
results, thus providing a consistency check that adding 
metallic iron can explain the infrared colours of other AGB 
stars as well. From the figure it becomes also clear that 
modifications of the iron mass fraction probably improves 
the result. To determine the iron mass fraction, detailed 
full radiative transfer calculations are necessary. 
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Table 3. Absolute magnitudes and infrared colour index 
as a function of mass-loss rate, for modelled SEDs with 
4% of the dust mass contained in non-spherical metallic 
iron grains. 



M (Mr, vr^M 




Ml 


L - fl2l 


4('— 81 


—8.8 


— 7.4 


1.4 


7('— 81 


—9.0 


—7.4 


1.6 


11-71 


—9.1 


—7.4 


1.7 


21 — 71 


—9.4 


—7.5 


2.0 


41 — 71 


—9.9 


—7.5 


2.4 


7(-7) 


-10.3 


-7.5 


2.8 


l(-6) 


-10.6 


-7.6 


3.0 


2(-6) 


-11.0 


-7.7 


3.3 


4(-6) 


-11.5 


-7.9 


3.6 


7(-6) 


-11.6 


-7.9 


3.7 


l(-5) 


-11.7 


-8.0 


3.7 


2(-5) 


-11.9 


-7.8 


4.1 


4(-5) 


-11.9 


-7.4 


4.4 


7(-5) 


-11.9 


-6.8 


5.1 


l(-4) 


-11.8 


-6.3 


5.5 


2(-4) 


-11.6 


-4.9 


6.7 


4(-4) 


-11.2 


-3.2 


8.0 


7(-4) 


-10.7 


-1.7 


9.0 


l(-3) 


-10.1 


-0.7 


9.4 



6. Discussion 

The results presented in this work provide the first sound 
identification of metallic iron in the circumstellar dust 
shell around AGB stars, without violating abundance con- 



straints. Demyk et al. (2000) have tried to include metallic 
iron in their model calculations of evolved stars, but be- 
cause they did not consider non-spherical grains, metallic 
iron could not be a bundant enough to account for the 
missing NIR opacity. Harwit et al. (2001 ) have concluded 
that metallic iron is necessary to fit the spectrum of VY 
CMa, although all Fe-atoms should be incorporated in 
metallic iron, implying that the amorphous silicate is com- 
pletely Mg-rich, and even then a large overabundance is 
required. 

From meteoritic and theoretical studies it is known 
that metallic iron can condense in a cooling gas of solar 



composition (e.g. Grossman & Larimer 1974; Lewis & Ney 
ozasa fc Hasegawa 1988[ ). In that case, metallic 
iron forms almost simultaneously with silicates, since its 
condensation temperature is only 50 - 100 K below that 
of silicates. In an oxygen-rich chemistry, metallic iron is 
stable above 700 K. Belo w that temp erature it will react 
to form FeS and/or FeO ( |jones 1990[ ). 

An alternative is that it will remain metallic because it 
is protected from the oxidizing environment by inclusion 
in a grain. This can be achieved through the formation 
of iron islands on the surface of silicate grains dGail fc| 



G 



10" 



10" 



^ 10" 



10 




10" 



4 6 
L-[12] 
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Fig. 6. The mass loss rate of oxygen-rich AGB stars as a 
function of the L — [12 /im] colour index, arising from ra- 
diative transfer calculations, indicated with the asterisks. 
The fit parameters of OH 127.8-1-0.0 are indicated with a 
square. The diamonds represent the objects from the sam- 



ple of Sylvester et al. (1999) with error bars on the mass- 
loss rate. The gas/dust ratio is assumed to be 100. With 
increasing mass-loss rates these are: o Get, WX Psc, CRL 
2199, OH 104.9-f2.4, OH 26.5+0.6 and OH 32.8-0.3. The 



mass-loss rate of o Get is taken from Loup et al. (1993). 



to the surface of the grain to form iron islands and stimu- 
late the condensation of additional Fe-atoms, to continue 
the growth of the islands. [Gail fc Sedlmayr (1999 ) sug- 
gest that the islands will eventually be covered by younger 
silicate layers, thus creating silicate grains with platelet 
shaped, i.e. non-spherical, metallic iron inclusions. These 
inclusions greatly increase the opacity of the grains in the 
NIR region. The model calculations presented in this work 
use separate grain populations, rather than metallic iron 
inclusions in silicate grains, but still provide an idea of the 
mass fraction contained in the iron inclusions. This result, 
together with the inclusion of crystalline water ice and the 
use of GDE for all grain shapes provide an important step 
in disentangling the dust composition of the circumstellar 
shell of AGB stars. No exotic assumptions were required 
to match the spectrum. A spherical outflow of material at 
constant velocity suffices to produce a good fit to the SED 
if chemical and shape properties of the grains are treated 
in sufficient detail. 

The presence and formation of metallic iron in astro- 
physical environments has been subject to many studies 



Sedlms yr 199S|). During the formation process of silicates 
Fe-atoms in the lattice will be replaced by the thermody- 
namically more favourable Mg. The Fe-atoms will migrate 



in the past (Jones 1990, and references herein). Studies 
of the interstellar extinction toward the galactic centre 
(GG) revealed a discrepancy in the 3-8 /im region with 
the standard dirty silicate opacity ( Lutz 1999| ; Lutz et al 



1996). Of course, this is an ISM line-of-sight, implying 
that amorphous carbon is an important dust component. 
On the other hand, the opacity contribution of amorphous 
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carbon is already included in the optical constants of the 
dirty silicate, i.e. a graphite-silicate mixture. Thus also in 
this case, non-spherical metallic iron grains might explain 
the missing opacity in the NIR region. In addition, it is 
known that superparamagnetic inclusions, such as metallic 
iron, in elongated silicate grains can cause these grains to 



align when a magnetic field is present (Mathis 1986). The 
alignment of silicate grains explains the polarization of 
starlight passing thro ugh the interstellar me dium, a sug- 
gestion first made by 3pitzer fc Tukey (1951 ). 



This theory appears to be supported by studies of in- 
terplanetary dust grains, which are collected by aircraft 
in the upper atmosphere of the Earth. Some of the col- 
lected particles are believed to be very pristine and prob- 
ably even from interstellar origin, rather than being re- 
processed during the formation of the solar system. These 
pristine particles consist of an amorphous matrix mate- 
rial with metal inclusions, and are referred to as GEMs 
(glasses with embedded metal and sulphides). Metallic 
iron inclusions ar e quite common in these sup posedly in- 
terstellar grains ( Bradley 1994 ; Martin 1995| ). Whether 
these grains really originate from the outflows of evolved 
stars remains subject to speculation at this moment. 
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